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0.   INTRODUCTION 

The  purpose  of  this  paper  is  to  present  and  analyze  a 
model  of  the  use  of  various  strategies  in  landing  patterns 
when  trying  to  recover  a  squadron  of  aircraft  returning  to  an 
aircraft  carrier.   Section  1  contains  a  brief  description  of 
the  landing  procedure  and  why  a  problem  arises.   Section  2 
describes  the  basic  model  from  which  we  derived  basic  equations 
of  interest.   In  Sections  3,  4  and  5  we  develop  expressions  for 
three  different  measures  of  effectiveness,  the  expected  time 
to  recover  all  aircraft,  the  variance  of  this  time,  and  the 
total  expected  fuel  used  in  recovery,  respectively.   The 
expressions  each  require  the  determination  of  certain  probabilities 
before  they  can  be  used  for  numerical  calculations.   Expressions 
for  calculating  these  probabilities  are  given  in  Section  6. 
Finally,  in  Section  7,  we  present  numerical  calculations  to 
demonstrate  how  the  model  can  be  used  to  investigate  various 
landing  strategies. 


1.   CARRIER  LANDING  PROCEDURES 

Figure  1  is  an  illustration  of  how  carrier  landings 
are  performed.   The  aircraft  (usually  about  30)  are  marshalled 
in  a  holding  pattern  20  nautical  miles  behind  the  carrier, 
separated  by  1000  feet  in  elevation,  starting  at  5,000  feet. 
Each  aircraft  is  given  a  precise  time  and  position  at  which  to 
leave  the  holding  pattern  and  follow  a  glide-path  to  the  landing 
deck  of  the  carrier.   These  times  and  positions  are  such  that 
nominally  the  aircraft  reach  the  landing  point  at  intervals  of 
1  min  15  sees.   Thus  with  30  aircraft,  if  we  count  time  from 
when  the  first  aircraft  reaches  the  landing  point,  and  if  all 
aircraft  land  successfully  on  their  first  attempt,  it  would  take 
36  min  15  sees  to  recover  all  the  aircraft. 

A  problem  arises  because  with  carrier  landings  not  every 
attempt  is  successful.   There  is  some  probability  that  an  air- 
craft is  either  waved  off  or  misses  the  arresting  cable.   We 
say  that  an  aircraft  which  does  not  land  successfully  on  a  given 
attempt  bolters.   Aircraft  which  bolter  join  a  queue  at  1,200  feet 
to  await  a  further  landing  attempt.   This  queue  we  call  the  bolter 
queue.   An  aircraft  in  the  bolter  queue  cannot  make  another 
landing  approach  until  there  is  room  for  it  in  the  landing  pattern 
If  all  aircraft  are  scheduled  so  that  they  arrive  at  the  landing 
point  every  1  min  15  sees,  then  aircraft  in  the  bolter  queue 
must  wait  until  all  aircraft  have  made  their  initial  attempt 
before  they  can  make  their  second  or  subsequent  attempt. 
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To  prevent  undue  delays  in  the  bolter  queue  it  is  common 
practice  to  leave  "bolter  holes"  in  the  landing  pattern.   A  hole 
is  simply  a  1  min  15  sees  time  in  which  there  is  no  aircraft 
scheduled  to  land  from  the  holding  pattern.   Figure  2  illustrates 
a  landing  strategy  with  eight  aircraft,  with  two  holes  left  in 
the  fourth  and  eighth  positions.   These  positions  could  be  filled 
by  aircraft  from  the  bolter  queue  if  necessary.   If  the  bolter 
queue  is  empty  when  a  hole  appears,  then  the  hole  remains  unfilled. 

The  probability  that  a  landing  attempt  is  successful 
depends  on  many  factors,  such  as  time  of  day  (day  or  night  land- 
ing) ,  weather,  moon  state  if  landing  at  night,  sea  state,  and 
experience  of  the  squadron  at  sea.   For  a  seasoned  squadron  which 

has  been  at  sea  for  some  time,  landing  with  calm  seas  on  a  clear 
afternoon,  the  probability  of  a  successful  landing  will  be  high. 
But  with  a  squadron  only  recently  at  sea  landing  in  rough  weather 
at  night  with  no  moon  or  in  rain,  the  probability  can  be  quite 
low. 

Before  formulating  the  model  a  few  further  points  of 
clarification  should  be  made.   Readers  may  think  (as  the  author 
originally  did)  that  the  problem  is  one  of  sequential  decision 
making.   If  an  aircraft  bolters,  then  one  should  update  all  the 
landing  times  of  the  aircraft  still  in  the  holding  pattern  in 
order  to  reschedule  the  boltered  aircraft.   This  may  be  an 
interesting  model,  but  it  cannot  be  implemented  in  practice  for 
many  reasons,  including  safety  and  an  unacceptable  increase 
in  radio  traffic.   Thus  the  problem  is  to  decide  a  ph.i-Oft.-L   what 
sequence  of  scheduled  landings  and  holes  to  use  for  a  given 
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A  -  SCHEDULED  AIRCRAFT  FROM  HOLDING  PATTERN 
H  -  BOLTER  HOLE 


H    A    A    A    H 


Time    (  1  unit  s  1  min,  15  sees  ) 
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Figure  2 :  Example  of  a   Landing  Strategy  with 
8   Aircraft. 


landing  probability  and  values  of  other  model  parameters.   The 
author  has  written  a  related  paper  (Marshall  [2])  which  deals 
with  the  model  building  problems  as  opposed  to  the  technical 
ones  of  the  final  model  used  which  are  described  below. 

Those  familiar  with  the  real  landing  problem  will 
recognize  some  simplifying  assumptions  which  have  been  made. 
It  is  assumed  in  this  paper  that  all  aircraft  eventually  land; 
none  are  sent  to  shore  or  to  aerial  tankers  for  refuelling. 
This  point  is  discussed  in  Section  7.   Also  it  is  assumed  that 
if  a  hole  is  present  and  the  bolter  queue  is  not  empty,  it  is 
always  possible  to  place  one  of  the  boltered  aircraft  in  the  hole 


and 


with 


n(t)  =  n(t-l)  +  q±  (1) 


n(2)(t)  =  n(2)(t-l)  +  2qin(t-l)  +  q±    ,  (2) 


n(0)  =  n(2)  (0)  =  qi 


Now  let   t   be  the  time  of  a  bolter  hole.   We  assume  that 
if  the  bolter  queue  is  not  empty  at   t-1,  an  attempt  is  made 
to  land  an  aircraft  from  this  queue  in  the  hole  at   t.   If  it 
does  not  land  it  returns  to  the  bolter  queue.   Thus 

N(t)  =  N(t-l)  -  U(t)  , 
where 

U(t)  =1   if  both   N(t-l)  >  0   and  the  aircraft  lands, 
=  0   otherwise. 

Note  that  for  this  case   t   cannot  be  zero.   Since 


P[U(t)  =  1]  =  p2  P[N(t-l)  >  0]  , 
we  have 


n(t)  =  n(t-l)  -  p2  P[N(t-l)  >  0]  ,  (3) 


n(2)(t)  =  n(2)(t-l)  -  2p2  n(t-l)  +  p2  P[N(t-l)  >  0]  , 


for  the  case  where   t   is  a  bolter  hole 


(4) 


Equations  (l)-(4)  are  used  in  the  next  three  sections 
to  determine  expressions  for  the  three  measures  of  effectiveness 
described  in  Section  1. 
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3.   .EXPECTED  RECOVERY  TIME  (r  ) 

Let   R  be  the  time  from  t  =  0   to  when  the  last  air- 
craft is  successfully  recovered.   For  a  given  strategy   S   with 
K  elements  and  with   m   aircraft  to  be  recovered,  the  smallest 
value  that   R  can  be  is   m  +  K  -  1.   Let   B   be  the  time  from 
t  =  m  +  K  -  1   until  the  bolter  queue  is  empty.   Then  clearly 

R  =  m  +  K-1  +  B,  (5) 

and  on  taking  expected  values,  if  we  let   r  =  E[R]   we  obtain 

r  =  m  +  K  -  1  +  E[B] .  (6) 

Since  it  will  take  a  geometric  number  of  attempts  to  recover 
each  aircraft  from  the  bolter  queue  after   t  =  m  +  K  -  1   with 
success  probability   p_,  we  have 


E[B]  =  n(m  +  K  -  1}  .  (7) 

P2 


For  a  given  strategy   S  =  {s,  ,  s_,  ...  ,  s„}   we  now  use 
equations  (1)  and  (3)  to  solve  for   n(m  +  K  -  1) .   First  define 


a.  =  P[N(s.  -  1)  >  0],        i  =  1,2,. ..,K  .       (8) 
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Thus   a.   is  the  probability  that  the  bolter  queue  is  non-empty 
when  the  i-th  hole  appears.   The  structure  of   a.   is  described 
later,  in  section  6,  but  notice  that  it  depends  on   s, ,  s»,...,s. 
and  on   p,   (and   p2   when   i  >  1) .   Using  (8) ,  (1)  and  (3)  it 
is  straight  forward  but  tedious  to  show  that 

K 
n(m+K-l)=mq,-p0   £  a.  .  (9) 

1    A   i=l  1 

Writing   r   as   r(S)   to  indicate  its  dependence  on  the  strategy 
S,    and  using  (9),  (7)  and  (6)  we  find  that 

m(p   +  q  )         K 

r(S)  = 1  +   J   (1  -  a.)  .  (10) 

P2  i=l 


In  order  to  use  (10)  for  calculating   r,  we  need  to  be 
able  to  determine  the   a.,  i  =  1,...,K.   Expressions  for  these 
are  derived  below  in  Section  6.   But  we  do  not  need  to  calculate 
them  to  determine  the  strategy  which  minimizes   r.   Since   a. 
is  a  probability,  all  terms  in  the  sum  on  the  right-hand  side 
of  (10)  are  non-negative,  and  this  sum  is  the  only  part  of  the 
expression  which  is  a  function  of  the  strategy  used.   Thus  to 
minimize   r  we  set   S  =  0.   That  is,  leave  no  bolter  holes  in 
the  landing  pattern. 

Equation  (10)  has  an  interpretation  which  makes  obvious 
this  result.   The  first  two  terms  in  (10)  give  the  total  expected 
number  of  attempts  it  will  take  to  land  m  aircraft  (the  -1 
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comes  from  our  choice  of  starting  at   t  =  0   with  a  scheduled 
landing) .   The  summation  term  gives  the  expected  number  of  empty 
holes  which  will  occur  using  strategy   S.   Since  there  is  always 
a  positive  probability  (often  very  small  in  practice  as  we  shall 
see)  that  a  bolter  hole  will  not  be  filled,  it  is  optimum  not 
to  leave  any. 
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4.   VARIANCE  OF  THE  RECOVERY  TIME  (v) 

An  important  factor  in  any  process  is  often  its  vari- 
ability.  Uncertainty  as  to  the  time  to  recover  all  aircraft 
causes  problems  in  the  scheduling  of  other  functions  aboard  the 
carrier.   Thus  it  is  important  that  we  know  how  the  variance  of 
the  recovery  time  changes  with  the  strategy   S.   It  may  be  possible 
to  choose  a  strategy  which  gives  significant  reduction  in  recovery 
time  variance  without  increasing  the  expected  recovery  time 
substantially. 

From  equation  (5),  if  we  define   v   to  be  the  variance 
of   R,  then 

v  =  Var[B]  .  (11) 

Thus  all  the  variance  in  the  recovery  time   R   is  in  the  recovery 
from  the  bolter  queue  after  the  scheduled  landing  pattern  has 
been  completed. 

A  simple  way  to  proceed  in  the  evaluation  of  v   in  (11) 
is  to  use  the  conditional  variance  relation 

Var[B]  =  E[Var[B|N(m  +  K  -  1)  ]  ]  +  Var[E(B|N(m  +  K  -  1)  ]  . 

Observing  again  that  each  aircraft  in  the  bolter  queue 
at  the  end  of  the  schedule  will  take  a  geometric  number  of  attempts 
to  land,  then 
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Var[B]  =  (n(m  +  K  -  l)q   +  Var[N(m  +  K  -  DD/p^  •       (12) 


An  expression  for  n(m  +  K  -  1)  is  given  in  (9) ,  but  we  need 
to  determine  one  for  Var [N (m  +  K  -  1) ] ,  or  equivalently ,  for 
n(2) (m  +  K  -  1) . 

Recall  that  equations  (2)  and  (4)  in  Section  2  give 

(2) 
recursive  expressions  for   n    (t)   when  there  is  a  scheduled 

landing  or  a  bolter  hole  at   t   respectively.   Straight  forward 

but  tedious  algebra  using  (2)  and  (4)  with  a  given  strategy 

leads  to  the  expression 


,  m+K-2  K                K 

n(  '   (m+K-1)  =mq,+2q  J    n(t)-2(q  +p9)  [  n(s.-l)+p    I    a...     (13) 

L       x  t=0          x  z  i=l    1      z   i=l  x 

By  using  (1)  and  (3)  we  obtain 


m+K_2  tm-1)mn  K  K 

I 
t=0 


(m-l)mq,        K 

I         n(t)  =  5 ^  +  q,  I     (s  +l-i)-p    I     (m+K-l-s.)a..   (1 

=0  i=l  i=l 


4) 


and 

K  K  K 

T   n(s.  -  1)  =  q-  I    (s.  +l-i)  -p,  I       (K-i)a.  . 
i=l     x        -1  i=l   1  Z    i=l 

From  (15),  (14),  (13),  (12)  and  (11),  after  some  algebra  one 
finds 


(15) 
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mq  (1  +  q  --  q  )    2q    K 
v(S)=— - ~-^ -±  I     (S.  +l-i)  (1-ct.  ) 

?2  Pl  i=1 


K  K 

+  2  J  (K-  i  +  |)a.  -  (  I   a.)2  ,  (16) 

i=l       z   1    ±=1      1 


where  we  have  written   v (S )   to  remind  the  reader  of   v's 
dependence  on  the  strategy  S.   Equation  (16)  can  be  used  to 
calculate  the  variance  of   R  once  the  probabilities   a., 
i  =  1,2,... ,K      have  been  determined  for  the  given  strategy. 

It  is  not  clear  from  (16)  how  the  variance  changes  with 
the  strategy.   We  note  that   v(0)   is  given  by  the  first  term 
which  is  simple  to  evaluate.   Note  also  that  if  we  replace   a. 

by  1.0  for  every   i,  the  second  term  drops  out  and  the  last 

2 
two  terms  are  each  equal  to   K   and  cancel.   Thus  for  strategies 

where  it  is  almost  certain  there  will  be  an  aircraft  in  the 

bolter  queue  when  a  hole  occurs,  there  is  virtually  no  reduction 

in  the  variance  over  the  no-hole  strategy. 

One  can  argue  that  a  strategy  can  be  found  that  makes   v 

arbitrarily  small  (at  the  expense  of  making   r   arbitrarily 

large) .   For  example,  suppose  we  use  the  (completely  unrealistic) 

strategy  of  leaving  20  holes  after  each  scheduled  landing. 

With  a  realistic  value  of   p,   the  probability  of  having  no 

aircraft  in  the  bolter  queue  at  the  end  of  the  schedule  would 

be  extremely  close  to  1.0.   Thus  we  would  have  virtually  no 
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variance  in   B,  and  hence  in   R.   The  point  of  this  example 
is  to  show  that  searching  for  a  strategy  which  minimizes  variance 
would  be  futile.   However  one  can  use  (16)  together  with  (10)  to 
investigate  tradeoffs  between   r   and   v   for  various  strategies. 
This  is  illustrated  computationally  in  Section  8. 
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5.   EXPECTED  FUEL  CONSUMED  IN  RECOVERY  (f) . 

An  important  factor  in  determining  which  landing  strategy 
to  use  is  that  aircraft  use  fuel  at  a  considerably  higher  rate 
in  the  bolter  queue  at  1200  feet  than  in  the  holding  pattern  at 
elevations  greater  than  5000  feet.   In  this  paper  we  assume  that 
<J)   is  the  ratio  of  the  fuel  consumption  per  unit  time  in  the 
bolter  queue  to  that  in  the  holding  pattern.   This  ratio  varies 
somewhat  with  aircraft  type  but  we  assume  that  it  can  be  esti- 
mated for  a  given  application.   Typical  values  for  modern  Navy 
aircraft  are  in  the  range  1.3-1.5.   Rather  than  work  with  units 
of  actual  fuel  consumption  (lbs/minute,  and  total  lbs  used  in 
the  recovery) ,  we  work  in  units  of  fuel  consumed  in  one  time  unit 
in  the  holding  pattern. 

From  equation  (5)  we  see  that  there  are  two  distinct 
parts  to  the  total  recovery  time,  the  scheduled  pattern  of  fixed 
length  m  +  K  -  1,  and  the  bolter  queue  recovery  time  following 
this  of  random  length   B.   We  calculate  the  total  expected  fuel 
consumed  by  considering  these  parts  separately.   Furthermore, 
we  break  down  the  fuel  used  in  the  fixed  interval   (0,  m  +  K  -  1) 
into  that  used  in  the  holding  pattern  and  that  used  in  the 
bolter  queue. 

Total  Fuel  Used  in  the  Holding  Pattern 

Since   t  =  0   is  the  time  that  the  first  aircraft  reaches 
the  landing  point,  this  aircraft  uses  no  fuel  in  the  holding 


18 


pattern  during  the  measured  time  interval.   The  second  uses 
1  unit  if   t  =  1   is  not  a  bolter  hole.   In  fact  for   0  <  i  <  s  , 
aircraft  (i  +  1)   uses   i   units.   in  general,  the  aircraft  which 
is  scheduled  to  land  at  time   t   uses   t   units  in  the  holding 
pattern.   If   h   is  the  total  fuel  used  in  the  holding  pattern, 
then  for  strategy   S   it  is  easily  shown  that 


h<s)  _  <m  +  K  -  !)(„  +  K)  _   5 

2  J   1     1 


K 

I 
i=l 


For  the  no-hole  strategy  we  have   h(0)  =  [m(m-l)]/2.   Using  this 
we  can  rewrite   h   in  the  form 


h(3)  .  &<»-i>. +  *(*» y-p  _  j  (17) 


K 

I 
i=l 


This  expression  is  independent  of  p,  because  once  a  strategy 
is  fixed,  the  actual  realization  of  successful  landings  has  no 
effect  on  h. 

Expected  Fuel  Consumed  in  the  Bolter  Queue  During  the  Landing  Schedule 

Because  of  our  choice  of  units,  the  total  expected  fuel 
consumed  in  the  bolter  queue  during  the  landing  schedule  of  length 
(m  +  K  -  1)   is  equal  to   <j>   times  the  total  expected  number  of 
aircraft-periods  spent  in  the  bolter  queue  during  this  period. 
Let   a   denote  this  total  expected  fuel.   Then 
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m+K-2 
a(S)  =  <J>    j;   n(t)  , 
t=0 


and  by  using  equation  (14)  we  obtain 


a(5)  = 


(m-l)mq,  <$> 


-*         K  K 

i~  +  q^      I       (S;L  +  1  -  i)  -  p2cj)  J    (m  +  K  -  1-  Si)  a±    . 


i=l 


i=l 


(18) 


Note  that  the  first  term  gives   a(#)  =  [ (m-ljmq,^ ]/2 . 

Expected  Fuel  Consumed  in  the  Bolter  Queue  Recovery  Time 

Recall  that  at  time   (m  +  K  -  1)    the  number  in  the 
bolter  queue  is  random  and  denoted  N(m  +  K  -  1) .   A  simple  way 
to  calculate  the  expected  fuel  consumed  in   B   is  by  using 
absorbing  Markov  Chain  results  with   N(m  +  K  -  1)   as  the  random 
starting  state  (see  Kemeny  and  Snell  [  l ]  ,  Chapter  3)  . 

Suppose  we  count  periods  starting  with  0  at  time  m  +  K  -  1, 
after  which  we  start  emptying  any  remaining  aircraft  from  the 
bolter  queue.   The  number  of  aircraft  in  the  bolter  queue  at 
each  time  after  this  forms  a  state  of  a  finite-state  Markov  Chain 
with  absorbing  state  0  and  transient  states  1,  2,     ...  ,  m.   The 
transient  part  of  the  one-step  transition  matrix  is  given  by 


Q  = 


(19) 


^2  J 
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which  is  an  m  x  m  matrix  with   q9's  on  the  main  diagonal,  p0 ' s 

on  the  lower  diagonal,  and  zeros  elsewhere.   If  we  start  in 

some  state,  say   n,  the  total  expected  fuel  consumed  will  be 

the  total  expected  number  of  visits  to  each  state  starting  in 

state   n,  multiplied  by  the  value  of  the  state.   The  expected 

number  of  visits  to  each  state  for  each  given  starting  state  is 

-1 


given  by  the  fundamental  matrix   (I  -  Q) 
that  for   Q   given  by  (19)  we  get 


It  is  easy  to  show 


(I  -  Q) 


-1 


P2 


1 
1 


1   1 


(20) 


an  m  x  m  matrix  with   l/p2   in  each  position  on  and  below  the 
main  diagonal,  and  zeros  above  the  main  diagonal. 

Now  let   C   be  the  total  fuel  consumed  in  emptying  the 
bolter  queue  with  expected  value   c.   Note  that  whenever  the 
chain  is  in  state  i,  there  are   i   aircraft  in  the  bolter  queue 
These  consume   <J>i   fuel  in  a  given  period.   If  we  write   N   for 
N(m  +  K  -  1)   to   simplify  the  notation,  then  from  (20) , 


E[C|N]  =  -£ 
P2 


V  (N2  -  N) 

I       (N  -  i)  =  —^r-^-   *• 


2p. 
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Using  equations  (9) ,  (13) ,  (14)  and  (15) ,  after  removing  the 
condition  on   N  we  obtain 


mq    [2    +     (m-l)q,  ]  <J>  K 

c(S)    =      — i j- -qi<|>    I     (s      +    1    -    i) 

2p2  1    ^       l 


K 

-  <j>    I    [q1(m-l-si+l)    -   p    (K-i)]    a±    .  (21) 

i=l 


Note  that  the  first  term  gives  the  expected  fuel  consumption  for 
the  no-hole  strategy,  c (0) . 

Let   f   be  the  total  expected  fuel  consumed  in  the 
recovery  of  all  m  aircraft.   If   f  (0)   represents  this  value 
for  the  no  hole  strategy,  then  from  (17) ,  (18)  and  (21)  we  get 

f  (0)  =  [m(m-l)  (p2  +  q1<(>(q1  +  P2)  )  +  2mql(j)]/2p2  , 

and  (22) 

K 
f  (S)  =  f  (0)  +   5!   (m+  i  -  s.  -  1)  (1  -  (|)(q  +  p  )a.)  . 
i=l         1  lax 

The  expressions  (22)  are  used  to  calculate  the  total 
expected  fuel  consumed  once  the   a.'s,  i  =  1,2, ...,K   have  been 
calculated  for  a  given  strategy.   Note  that  if   p,  =  p_   and 
4>  =  1.0,  then  all  terms  in  the  sum  in  (22)  are  non -negative . 
Thus  if  the  fuel  consumption  in  the  bolter  queue  were  the  same 
as  that  in  the  holding  pattern,  and  if  the  probability  of  landing 
were  the  same  for  each  queue,  then  the  optimal  strategy   would 
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be  to  leave  no  holes  when  trying  to  minimize  total  expected 
fuel  consumed. 

In  Section  7  equations  (10) ,  (16)  and  (22)  are  used 
in  numerical  computations  to  investigate  various  landing 
strategies.   All  these  require  computation  of  the   a.  's,  and 
these  are  investigated  in  the  next  section. 
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6.   PROBABILITY  CALCULATIONS 

Before  numerical  calculations  can  be  carried  out  with 
the  expressions  derived  so  far,  we  need  to  be  able  to  evaluate 
the   a.'s,  i  =  1,2,...,K.   Of  course  for  the  no-hole  strategy 
there  are  none  to  calculate,  so  we  assume  that   K  >  1.   In  this 
section  we  present  an  algorithm  to  evaluate  the  probabilities 
recursively. 

Consider  a  strategy   S  =  {s..,  s~  /  ...  ,    sK)   for  a 
given   K  >  1,  where   s,   is  the  kth  hole.   Now  define 


&±   =  s1  -  1, 

d1=s,-s11-l,       2<k<K, 
k    k    k-1     '         —   —    ' 


D.  =  I      d.  =  s,  -  k,     l<k<K. 
K    i=1   l     K  -    - 


Also  define  the  binomial  probability  mass  function 

b(j;n,p)  =  (*)  pj(l-p)n"j,    0  <  j  <  n 
=  0  n  <  j  . 

Finally  define 

g.  (j)  =  P[N(s,-l)  =  j],       0    j    D, 


and 


hR(j)  =  P[N(sk)  =  j]  ,        0  <  j  <  Dk 
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Thus   9i,(J)   gives  the  probability  of  having   j   aircraft  in 
the  bolter  queue  immediately  preceding  the  k-th  hole,  and 
h,  (j)   gives  a  similar  probability  immediately  after  the  k-th 
hole.   Clearly   aR  =  1  -  gk(0),  1  £  k  <_  K. 

The  following  algorithm  can  be  used  to  calculate  the 
g,  •  s  recursively: 

Algorithm 

1.   Set   gi(j)  =  b(j;d1,q1) ,  0  £  j  £  D± 

and      k  =  1 . 


2.   Set 


Mj)  = 


gk(j)  +  gk(j  +  Dp2  /  j  =  o  , 

gk(j)q2  +  gk(j  +  Dp2  ,         i  £  j  £  Dk_1   , 
gk(j)q2  /  j  =  Dk  > 

0  Dk  <  j  . 


3 .   Set 


gk+i(^  =  X  \U-u  b(i;dk+i^i}   '         °  £  J  £  d 

1=1 


4.  If   k  +  1  <  K,  set   k   equal  to   k+1   and  return  to  Step  2. 
Otherwise  go  to  Step  5. 

5.  Set   a=l-gk(0),   1  <_  k  <_  K . 

6.  Stop. 

The  calculations  in  this  algorithm  follow  directly  from  the 
definitions  of  the  probabilities. 
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7.   NUMERICAL  EXAMPLES 

In  this  section  we  use  equations  (10),  (16)  and  (22) 
to  demonstrate  the  effects  of  different  landing  strategies  on 
the  expected  recovery  time   r,  the  variance  of  this  time   v, 
and  the  expected  fuel  consumed  in  recovery   f,  respectively. 
The  following  (typical)  parameter  values  are  used  in  all  cal- 
culations : 

Number  of  aircraft,  m  =  30, 
Bolter  queue  fuel  ratio,  <f>  =   1.4. 

We  also  restrict  ourselves  to  strategies  with  "equally  spaced 
holes."   For  example,  if  we  place  a  hole  at  every  eighth  position, 
we  would  have   S   =  {7,15,23,31},  and  from  the  definitions  in 
Section  6,  d.  =  7,  i  =  1,2, ...,K  with   K  =  4.   These  strategies 
are  ones  which  are  usually  used  in  practice,  and  the  appropriate 
choice  of   d.,  as   p,   and  p~   vary,  is  discussed.   The  value 

of   d.   gives  the  number  of  aircraft  between  holes   i-1   and   i. 

l   3 

Note  that  the  number  of  aircraft  after  the  last  hole  cannot 
exceed   d.,  and  is  equal  to  it  only  when   m   is  divisible  by   d.. 

In  the  first  set  of  calculations  we  assume  that  the 
probability  of  landing  successfully  from  the  bolter  queue  is 
the  same  as  that  from  the  holding  pattern; i .e .  p„  =  p, .   Table  1 
shows  the  values  of   r,  v,  and   f   for   d,  =  1 ,2 , . .  .  , 10 ,15 , 
and  the  "no-hole"  case   (d,  =  30),  for  values  of   p,  =  p_  = 
0.5,  0.6,  0.7,  0.8,  0.9.    The  results  are  plotted  in  Figures  3, 
4  and  5. 
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The  results  for   r   illustrate  our  previously  observed 
result,  that  it  is  minimized  by  using  the  no-hole  strategy. 
However,  one  can  see  from  the  curves  in  Figure  3  how  insensitive 
the  expected  recovery  time  is  to  the  strategy  used.   For  example, 
p.  =  p„  =  0.7   there  is  virtually  no  change  in  expected  recovery 
time  for  all  strategies  with   d.  >  4.   Notice  that  in  all  cases 
shown  there  is  a  sharp  increase  in   r   for  very  small  values  of   d. 

The  results  for   v   are  shown  in  Figure  4.   Again  we 

see  how  insensitive   v   is  to   d.   for   d.   larqe,  but  how 

sensitive  it  is  for   d.   very  small.   These  curves  illustrate 

the  point  made  at  the  end  of  Section  4 ,  that   v  can  be  reduced 

at  the  expense  of  increasing   r.   For  example,  if   p   =  p   =  0.6, 

the  no-hole  strategy  gives  an  expected  recovery  time  of  49.0 

units  with  a  variance  of  33.3.   With   d.  =  5   (a  hole  every  sixth 

i 

position),  r   increases  to  49.1  and  the  variance  is  reduced  to 
32.6;  both  are  insignificant  changes.   But  if   d.  =  1   (a  hole 
every  second  position),  then   r   increases  to  59.9  while   v 
is  reduced  to  only  5.4. 

The  results  for   f   are  shown  in  Figure  5 .   Note  that 
for   p.,  <  0.8   there  is  a  minimum  point  on  each  curve  (the 
nimum  for   p.,  =0.9   lies  outside  the  plotted  range  of 


mi 


d.).   These  minima  have  been  marked  for  clarity.   For 

i 

example,  if   p   =  p_  =  0 . 7   the  strategy  which  minimizes  the 
expected  fuel  consumed  is  to  leave  a  hole  in  every  fifth  position 
(d.  =  4),  or  S   =    {4,9,14,19,24,29,34}.   Note  that  with  this 
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strategy  there  is  very  little  increase  in  the  expected  recovery 
time,  and  a  small  decrease  in  the  variance,  when  compared  to 
the  no-hole  strategy.   Note  however  that  the  expected  savings 
in  fuel  is  only  23  units  (714-691),  or  3.2%  of  that  used  with 
the  no-hole  strategy.   Thus  the  optimum  "equally  spaced  hole" 
strategy  gives  little  improvement  over  the  simple  no-hole  strategy 
for   p,  =  p~  =  0 .  7,  m  =  30 ,  and   cj)  =  1.4.   The  reader  can  see 
from  Figure  5  that  this  result  holds  quite  generally. 

The  small  fuel  savings  seen  in  the  calculations  given 
in  Table  1  make  it  questionable  as  to  whether  or  not  it  is  worth 
searching  for  an  optimal  strategy  among  those  that  do  not  have 
equally-spaced  holes.   Although  this  could  certainly  be  done  using 
(22)  we  do  not  pursue  this  route,  but  consider  another  alternative 
not  yet  discussed  in  the  paper.   A  number  of  aircraft  are  now 
fitted  with  automatic  landing  systems,  which  when  operated,  make 
landing  virtually  certain.   There  are  a  number  of  reasons  why 
these  devices  are  not  always  used.   Constant  use  would  lead  to 
a  reduction  in  landing  skills  by  the  pilot.   This  would  lead  to 
a  dangerous  condition  in  case  of  failure  of  the  automatic  system, 
and  would  also  lead  to  pilot  discontentment  and  motivation 
problems.   A  reasonable  question  to  ask  is  what  strategy  should 
be  used  with  regard  to  the  automatic  landing  system. 

Equations  (10),  (16)  and  (22)  can  be  used  to  demonstrate 
the  effect  of  using  the  automatic  landing  system  when  making 
an  attempt  from  the  bolter  queue.   To  do  this  we  repeat  our 
previous  calculations  with  the  condition   p   =  p    replaced  by 
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Figure  3:  Expected  Recovery  Time  with  p  =  p  ,  m  =  30 
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Figure  4:  Variance  of  Recovery  Times  with  p  =  p  ,  m  =30 
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Figure  6:  Expected  Recovery  Time  with  p  =1.0,  m  =  30. 
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Figure  7 


Expected  Fuel  Consumed,  p  =  1.0,  m  =  30,  <f>  =  1.4 
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p   =  1.0.   The  results  are  given  in  Table  2,  and  are  plotted 
in  Figures  6  and  7. 

Figure  6  shows  the  expected  recovery  times  for  various 
values  of   p, ,  as  a  function  of  the  bolter  hole  spacing.   The 
same  general  shape  is  observed  as  was  found  in  Figure  3  for 
Py    =  p,  .   However,  for   p-.   small   (0.5  -  0.6)   there  is  a  con- 
siderable reduction  in  recovery  time  as  one  might  expect.   These 
curves  can  be  used  to  measure  this  reduction.   For  example,  if 
p,  =0.6   (possibly  from  bad  weather  conditions  at  night)  and  a 
strategy  is  used  with   d.  =5,  the  expected  recovery  time  can 
be  reduced  from  49.1  time  units  to  41.1,  a  reduction  of  10  minutes 
(16.3%)  by  using  the  automatic  landing  device  on  attempts  from 
the  bolter  queue. 

A  look  at  Table  2  will  show  the  reader  that  considerable 
reductions  in  the  variance  of  recovery  time  are  achieved  with 
p~  =  1.0.   Because  the  variances  are  so  small  in  all  cases,  they 
have  not  been  plotted. 

Figure  7  shows  the  expected  fuel  consumed  in  recovery 
for  various  values  of   p, ,  and  these  curves  should  be  compared 
with  those  in  Figure  5.   Notice  that  they  have  the  same  basic  shape 
and  that  the  minimum  points  have  been  marked.   The  minima  are 
more  pronounced  in  Figure  7,  showing  a  greater  percentage  savings 
in  fuel  over  the  no-hole  strategy.   For  the  case  described  above 
with   d.  =  4,  p   =  0.7,  m=  30,  (J>  =  1.4,  but  with   p   =1.0,  the 
expected  savings  in  fuel  is  54  units  (685-631),  or  7.9%  compared 


with  only  3.2%  for   p   =  p 
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Not  all  possible  strategies  have  been  analyzed  in  this 
paper.   We  have  presented  enough  to  illustrate  the  output  which 
can  be  obtained  from  the  model.   It  should  be  quite  straight- 
forward to  program  on  a  hand-held  calculator  the  relevant 
equations,  so  that  instant  calculations  can  be  made  aboard  a 
carrier  with  very  little  resources.   In  this  way  any  conditions 
prevalent  on  a  given  recovery  can  be  taken  into  account. 
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